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Abstract 

Taking a hint from Dirac's large number hypothesis, we note the existence of cosmic combined 

o 

CN \ conservation laws that work to cosmologically long time. We thus modify or generalize Einstein's 

theory of general relativity with fixed gravitation constant G to a theory for varying G, which can 

hj ; 

^ — , . be applied to cosmology without inconsistency, where a tensor arising from the variation of G takes 

<n : 

the place of the cosmological constant term. We then develop on this basis a systematic theory 
■ of evolving natural constants m e ,m p ,e,h,kB by finding out their cosmic combined counterparts 

in 

involving factors of appropriate powers of G that remain truly constant to cosmologically long 



time. As G varies so little in recent centuries, so we take these natural constants to be constant. 
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I. GENERAL RELATIVITY AND LARGE NUMBER HYPOTHESIS 



In Einstein's theory of general relativity as in Newton's theory of gravitation, the strength 
of the gravitational interaction is described by a fixed dimensional constant, namely the 
Newtonian gravitation constant Gn = 6.7 x 10 _11 m 3 kg _1 s _2 . Together with the velocity of 
light in vacuo c = 3 x 10 8 m/ s which serves to convert time into length by putting c = 1, there 
occur also in Einstein's theory G N /c 2 = 7.4 x 10~ 28 mkg _1 and G N /c 4 = 8.2 x 10~ 45 mJ -1 
which respectively serve to convert mass and energy into length by putting also Gn = 1- 
Einstein's theory has been applied to cosmology; in addition to the dimensional constants 
Gtv another cosmo logical constant with dimension of length to the power minus two, was 
once specially introduced by Einstein in his attempt to construct a static model for the 
universe. This soon has lost his favor as observations of the red shifts of extra-galactic 
nebulae show defintely that the universe is expanding, but as an additional parameter, 
cosmological constant is still used today in trying to fit the observational data with the 
non-static homogeneous cosmological model using the Roberston- Walker metric. Rough 
estimates have thus been obtained for the age t of universe at present to be ~ 10 9 years or 
Pd 10 10 years more recently. 

On comparing the ratio e 2 / (Gm p m e ) of the electrostatic to the gravitational force between 
the proton and the electron in an hydrogen atom, which is a large dimensionless number of 
the order of 10 39 , with the ratio t/(e 2 /m e c 3 ) of the age t of universe at present to the time 
needed for light to travel a distance of the classical radius of the electron, which is also a 
large dimensionless number of the order of 10 39 or 10 40 , one proposes .1] that a large number 
equation 

e 2 /(Gm p m e ) « t/(e 2 /m e c 3 ) (1.1) 

should hold. Further from Hubble's estimate p = (1.3 to 1.6) x 10 -30 g/cm 3 for the density of 
matter due to the extra-galactic nebulae averaged over cosmic space and a factor thousand 
or hundred times to include dust or invisible matter, one can estimate the ratio Mjm v of 
the total mass M of matter in the universe at present of radius R to the proton mass, which 
is roughly of the order 10 78 or 10 80 , one proposes |lj that another large number equation 

M/m p « [t/(e 2 /m e c 3 )] 2 (1.2) 
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should also hold. By a large number equation we mean wherein some unknown factors of 
small numbers are understood to be there but not written out. Eddington takes seriously 
the relation 

M/m H = N w [e 2 /(Gm p m e )] 2 w 10 78 (1.3) 

and advocates in his book "Fundamental Theory" that in the uranoid (model of universe 
with iV hydrogen atoms) when one studies one particular hydrogen atom the other N — 1 
hydrogen atoms must also be considered as representing the rest of the universe. 

Based on these empirical large number equations (1.1) and (1.2) Dirac proposes the 
large number hypothesis that during the evolution of the universe the strength G of 
gravitational interaction gradually weakens as time t goes on according to G oc t^ 1 while 
mass is continuously being created and so increasing with time according to M oc t 2 . He 
reached at this hypothesis by adding an additional assumption that the atomic constants 
remain unchanged and do not evolve. One can, however, arrive at Dirac's large number 
hypothesis by the following evolving atomic constants, say, both m p and m e oc t 2 while 
e 2 a t 3 , leaving c alone not evolving. This corresponds to Eddington's assumption that N 
remains unchanged. 

In this paper we shall take the large number equations, though empirical at the present 
stage of knowledge, rather seriously. On dividing (1.2) by (1.1) we obtain a relation free 
from all atomic constants e,m e ,m p , 

GM ss ch (1.4) 

which exhibits clearly the inconsistency of applying Einstein's theory of general relativity to 
cosmology. If one regard this result, though empirical, as true, then clearly G and M cannot 
both conserve in cosmologically long time, and this is in contradiction to the assumptions 
inherent in Einstein's theory of general relativity that G be constant and matter do conserve. 
Hence one needs to modify Einstein's theory for application to cosmology. But we know, 
for phenomena including the crucial tests that occur, from the cosmic point of view, in a 
small neighborhood and for short duration, Einstein's theory is a good approximation. This 
has been confirmed by observations and experiments and must not be spoiled. We note, 
according to Dirac's large number hypothesis, it is neither M nor G but the combination 
G 2 M that conserves in cosmologically long time. We shall take such cosmic combined (c.c.) 
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conservation law working to cosmologically long time as a hint in our attempts to modify or 
rather to generalize Einstein's theory of general relativity, to be consistent with Dirac's large 
number hypothesis, and develop systematically a theory of evolving natural constants. In 
so doing we find that it is possible to generalize Dirac's large number hypothesis somewhat 
to G oc t~ n with n not necessarily equal to one. It seems that c. c. conservation laws can be 
more simply found and stated, but the variation of G, together with that of g^ u , as functions 
of space-time must be left to the solution of various peoblems. 

II. THEORY WITH VARYING G 

The problem of non-static homogeneous matter-dominated cosmological model is a simple 
case to work with mathematically. Because there is only one independent variable t, it is 
easier here to find out how natural constants shall evolve and how the usual conservation 
laws shall be combined to work in cosmologically long time. We shall not make Dirac's 
assumption that atomic constants do not change in the long. (In this paper by long or 
short we mean cosmologically long or short time and similarly by large or small we mean 
cosmologically large or small distance). But we do follow the spirit of Dirac's large number 
hypothesis and generalize it into 

G cc t~ n , (0 < n < 2), so t oc G~ 1/n (2.1) 

where n can be, but not necessarily, equal to one. Then (1.3) leads to 

G 1+1/n M oc t° i.e. G 1+1/n M = M conserves. (2.2) 

To save one's worry about the dimension of the cosmic combined (cc.) quantity like M, it 
would be most simple to introduce a constant but arbitrary dimensional conversion factor 
G to convert mass into length by putting it to be G = 1, so that in (2.2) and similar 
equations in the following we understand that the c. c. quantity (like M ) is always of the 
same dimension as the original quantity (like M). Thus the equation (2.2) is meant to be, 
when writen in full, 

(G/G ) 1+1/n M = M (2.2a) 

As usual we take Gn for Gq. In order to conform to such a c. c. conservation law for M we 
modify the corresponding action integral for simple matter without internal stress which in 
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Einstein's theory is given, as shown in Dirac's book |2( on general theory of relativity, by 

I m = - j{g ia ,^p v ) 1,2 d^x (E2.3a) with constraint f/\ M = (E2.3b) 
into a similar expression 

I m — — J (g^ u p fl p u ) 1 / 2 d 4 x (2.3a) with modified constraint p fJj , IM = (2.3b) 
Here with a change of notation 

= p»/(g al 3p a pP) 1/2 (E2.4a), so = 1 
(^pV) 1/2 = PmV(-9) so fpf = u» PmV /(-g) (E2.4b) 

the constraint in Einstein's theory gives the conservation of mass and may be written as the 
equation of continuity 

/^W = M (E2.4C) = (E2.4d) 

Similarly with a change of notation with the cosmic combined quantities 

(g af} p a pP) 1/2 = PmV(-9) (2-4a) = u"p my /(-g) (2.4b) 

where 

p m = G W/n pm (2 .5) 

the modified constraint (2.3b) in the modified theory gives the conservation of the cosmic 
combined mass 

J -M-O^M (2.4c)(2.5c) 
for homogeneous G = G(t) and the corresponding equation of continuity 

{PmU")-„ = (2.4d) 

In case of a mass point described with a three dimensional delta function factor in (E2.4c) 
and (2.4c), we obtain from (2.5c) the cosmic combined constant mass 

m = G 1+1/n m (2.5a) 

which holds for a proton, an electron, an atom or a molecue, contrary to Dirac's assumtion 
that atomic constants like m e and m p do not evolve. From this we see that our modification 



or generalization of Einstein's theory consists of multiplying the integrand of the source 
action integral by the dimensionless varying factor Q 1+1 / n and introducing new variables 
suitable for expressing the modified cosmic combined conservation law. This rule will be 
kept throughout for other action integrals involved in the comprehensive action principle, 
cf. §6, §8 and §12 later. 

As to Einstein's original gravitational action integral where G^ 1 is constant 



This is the most general expression for an action integral involving only field variables 
G, g^v without any dimensional constants that will reduce to Einstein's I g in case G be a 
constant. As the numerical factor (167r) _1 in Einstein's theory is determined by comparison 
with Newton's theory in the so called non-relativistic approximation (which is good for 
phenomena with velocity small compared with that of light), so the factor (167r) _1 in our 
modified theory is determined by comparison with Einstein's theory in the non-cosmologic 
approximation (which is good for phenomena of duration and distance short compared with 
that of the universe, as will be shown in §5). The numerical constant w will be determined 
from the non-static homogeneous matter-dominated cosmological model in §4 by using our 
generalized large number hypothesis (2.1). We prefer the choice of w = 8 for the case k = 
of the Robertson- Walker metric. 

In later sections we shall see that such a modification or generalization can be carried 
through other source action integrals like that for Maxwell's theory of electromagnetism 
(§6) or that for Dirac's theory of electron (§8). In §12 we consider statistical mechanics. 
We thus arrive at a systematic theory of evolving natural constants passing all the crucial 
tests of Einstein's theory of general relativity and agreeing with the empirical large number 
equalities (1.1), (1-2) or (1.4) when applied to cosmology. We leave the exact numerical 
value n introduced in the generalized Dirac's large number hypothesis to be determined in 
the future by fitting accurate relevant astrophysical or cosmological observations. For the 
moment one may take n — 1, the only integer in tne domain < n < 2, if one likes, 
following Dirac's original large number hypothesis. 




G~ 1 R° y /(-g)d 4 x where R* = R^ (E2.6) 



we multiply the integrand by the factor G 1+1 ^ n and add a kinetic term thus: 
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III. VARIATIONAL EQUATIONS 



It is convenient (cf.(2.6) and (2.1)) to introduce the dimensionless new variable 

= G 1/(2n) , (t) 2 = G l/n (3.1) 

so that (2.6) becomes 

I 9 = (167T)- 1 J (tfR^ - wh^g^y/i-g^x (3.2) 
We obtain the variation as 

6I g = (IGtt)- 1 J [-N^5 9lxv + 2<S>64>W(-g)d*x (3.3) 

with 

= <p 2 \Ri - + (<f>X - i<t> 2 )?X 

- W [0 ;a ;/3 -(l/2)0 ;(T 0^] (3.4) 

(here we use the Palatini identity to express SR^ in terms of the second covariant derivatives 
of Sg a /3 and integrate by parts twice), and 

^ = Ry + w^ (3.5) 

We note the identities (which can be shown from (3.3) by using an infmitesmal trans- 
formation of coordinates or explicitly verified from (3.4) and (3.5) by using the formula 
(A x ).^. u — (A x ). v .^ = — A a R x aiiv with \ = u summed), namely 



To calculate the variation of I m , (2.3a), we follow closely Dirac's treatment for I m , (E2.3a), 
as given in his book [2j on general theory of relativity. (Only we use the letter u instead of 
Dirac's v to denote the four-dimensional velocity) 

SI m = -(1/2) J f^6g^(-g)d 4 x - J Ufl 6p^d 4 x (3.7) 

with 

= {g a pp a P P )- 1/2 p"p V /^{-9) = P m uV (3.8) 
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by a change of notation that 

«" = (<m£V) _1/ V (3.9) 

PmV(-g) = (g«pp a p") 1/2 (3.io) 

The second part of 8I m , after substituting 

5p» = (p%»-p%n^ (3.H) so 5^=0 (3.12) 

in compliance with the above constraint (2.3b), where b^ is the four-dimensional virtual 
displacement of the element of matter, becomes after integrating by parts 

-J u^Vd A x = - J -tr\r), v #x 

= J(u^-u u ^)p u b»d A x (3.13) 

The vanishing of this integral for an arbitrary gives the geodesic equation of motion for 
the element of matter, namely 

p m u v {u^ v ) = p m u v {u^. u - u v .^) = pmU^U^ = (3.14) 

the other part vanishing identically because we have from (3.9) the identity u"u v = 1 so 
u v u v;il = 0. Here we note the identies similar to (3.6) 

= («\) ; , = {p m u»)^u v + PmU^U^ = (3.15) 

The identities (3.6) and (3.15) show that the system of variational equations obtained from 
the comprehensive action principle 5I to t = SI g + 5I m = 0, namely the gravitational field 
equations for g^ v and (/>, 

NZ = <t> 2 [K - imK€\ + (<pX - - - ^m«r€\ 

= -8vrTf = -8np m u a uP (3.16) 

$ = (^ + ^0) = O (3.17) 

and the equations of motion (3.14) for the elements of matter are compatiable but indeter- 
minate. This is only natural for invariant action integrals because the system of equations 
must allow g^ v to change with transformation of coordinates. 
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We note the following combination of (3.16) and (3.17) gives a particular simple equation, 

AT; + 0$ = ( w /2 - 3)(<j)% = -8vrf; = -8np m (3.18) 

It is interesting to write the field equations for in a form comparable with that in 
Einstein's theory involving a cosmological term. Substituting (3.4) into (3.16) and dividing 
both sides by <fi 2 we note that on the right-hand-side we have from (2.5) and (3.1) 

p m /<J) 2 = G Pm (3.19) or = GT% (3.20) 

Hence we obtain 

i^-(l/2)^ + A2 = -87rGr; (3.21) 
where on the right G is not a constant but varies as t~ n while on the left 



A£ = U 2 Y£ 



(.2\;er rfll / ,2 



- (l/2)0 ;(T ;CT ^]/0 2 (3.22) 

differs in its tensor character from Einstein's A5f = — Xdg and is originated from the 
variation of G , remembering (3.1). In case G = const in a region, vanishes in the same 
region. (Here we use the capital Greek A^ in (3.21) to compare with the similar equation in 
Tolman's book j^J on relativity, thermodynamics and cosmology which we write here with 
Gm restored and one index lowered 

-8irG N TP = Ri- (1/2)RX + A5f 
The corresponding equation in Weinberg's book 0] on gravitation and cosmology is 

Ri - (l/2)i£<$£ - A5f = -8ttG n TP 
so A = —A. To compare with the latter we should use = — A^) 



IV. MATTER DOMINATED COSMOLOGICAL MODEL 

We can adopt as usual the Robertson- Walker metric because in simplifying the metric 
to this form only considerations on symmetry and freedom of coordinate transformation 
have been used, but no use is made of the field equations. With t denoting the cosmic 
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time and r, 9, (p the dimensionless co-moving coordinates, ds 2 = dt 2 — R 2 (t){dr 2 / (1 — At 2 ) + 
r 2 d9 2 + r 2 sin 2 9d(f 2 }, k — 1, 0, or — 1 all the equations of motions for the element of matter 
are trivially satisfied by u A = w 4 = 1, u l = ui = i.e. by being at rest in the co-moving 
coordinates. From the modified constraint (2.3b) we obtain, using dots for time derivatives, 

{p m R 3 [(l - kr 2 )- 1 ^ 2 sin 9]Y = i.e. p m = p m (t )R 3 (t )/R 3 (4.1) 

where the constant to is really arbitrary but we often choose it to be the time or the age of 
universe at present. Then (3.18) can be integrated to give 

i? 3 (0 2 )' = (3 - w/2)- l ^p m {t )R\t Q )t (4.2) 

Here the constant of integration additive to t is chosen to be zero so that t is counted since 
the big bang when R = 0. We know according to the generalized large number hypothesis 

(2.1) that by (3.1) 2 oc t^ 1 in non-static homogeous cosmological model, hence (4.2) shows 
that R ext. Writing proportions as equalities 

2 = to0 2 (to)A (4.3) 
R/R(t ) = t/t (4.4) or R = (3 ex t (4.5) 

f3 ex being the dimensionless proportional constant of the expanding universe, we see that 

(4.2) is satisfied and gives 

87rp m (t o )/0 2 (to) = 87iG(t )p m (t ) = (w/2 - 3)/t 2 (4.6) 

Substituting (4.3) and (4.4a,b) into (3.17), (3.16), we see all the field equations are satisfied 
if we choose for our model with 

w = 8(1 + k/p 2 ex ) (4.7) 

For example the field equation (3.17), namely 

.2 .... 
R£(l> + w(f>% = 6[R/R+(R +k)/R 2 ](f) + w[<f> + 3(R/R)<j>] = (4.8) 

and the field equations iV* = N 2 — N% — which all coincide into, cf.(3.16), 

Nl = <P 2 [-(R 2 + k)/R 2 - 2R/R] - (0 2 )" - 2(<P 2 ) R/R 

. 2 

+w<f) /2 = (4.9) 
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both agree to determine w as above, while the equation 

. 2 

jV 4 4 = <f) 2 [-3(R +k)/R 2 ]-3((j) 2 yR/R 
-wifr /2 = -8npg (4.10) 

gives a relation similar to (4.6) but with the factor (w/2 — 3) replaced by (w/8 + 3k//3 2 x ) 
. These two factors, however, agree to be equal to (1 + Ak/(3 2 x ) by (4.7). Thus (4.3) and 

(4.5) are indeed the solution of our non-static homogeneous cosmologic model, where (see 

(4.6) ,(4.7)) 

8nG(t )p m (t )t 2 = 1 + 4k/ Pi = 8irGp m t 2 , (4.11) 
holds for arbitrary t or t. Together with 4irp rn R 3 /3 = M, we obtain from (4.11) 

GM = ((3 3 ex /6)(l + 4k/(3 2 ex )cH (4.12) 

in agreement with the large number relation (1.4) from where we started. With (2.1) and 

(2.2) , (4.12) becomes simply the conservation of M. Also (4.11) shows by (2.1) that p m cc 
£-(2-n) which suggests the domain for n given in (2.1) above. We have from (3.22) and 

(4.3) (4.4) , or cf. (4.9) and (4.10), 

A\ = A 2 , = Af = (1 + k/(3 2 ex )/t 2 
At = (2-k/p 2 ex )/t 2 (4.13) 

We believe that the action integral I g should be given once for all, being independent of 
the problems as here treated. From (4.6) we must choose w > 6, and this excludes by 

(4.7) the case of k = —1 for j3 2 x < 1. For the same reason the case of k — +1 can occur 
only for w > 16. Our first preference is to choose w — 8, thereby the case of flat three 
dimensional space k = is determined. The second choice w = 16, which corresponds to 
the case k — +1 and (3 ex — 1, is far less probable, we think. This question may be settled 
by noting the difference in (4.13), i.e. 

A4/AJ = 2 for w = 8 (4.13a) 
A 4 JA{ = 1/2 for w = lQ (4.13b) 

Of course a fresh analysis of the data of modern observation according to the formulae of 
the present theory would be very important. Our cosmology term varies inversely with the 
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square of the cosmic time. In terms of the Hubble function obtained from (4.4) (please note 
the different numerical coefficient from Einstein's theory) 

H{t) = R(t)/R(t) = 1/t (4.14) 

both the age of universe t and the values of our cosmological term at present can be 
determined from the Hubble constant H = H(t ) defined by small cosmic redshifts as 
discussed in §9. This seems to give a natural explaination for the order of magnitude of the 
cosmological constant at present. 

V. NON-COSMOLOGICAL APPROXIMATION 

We are intrerested in examine whether our modified theory may pass the crucial tests 
that support Einstein's theory of general relativity. Here we shall consider in particular that 
refering to the advance of perihelion of Mercury. We need only the exterior solution for 
g^ v and <fi good for, cosmologically speakng, a small region of space around the sun and a 
short inteval At around t say the present epoch. The motion of the planet e.g. Mercury 
will then be obtained from the geodesic equations of motion (3.4). We consider at first the 
cosmic background, i.e. the cosmological model studied in §4. We write the dimensionless 
co-moving coordinate in the Robert son- Walker metric in this section as tr-w and introduce 
a new variable r of dimension length 

r = R{t )r R _ w (5.1) 

The metric becomes for example in the case of k — 

ds 2 = dt 2 - [R 2 (t)/R 2 (t )}(dr 2 + r 2 d6 2 + r 2 sin 2 6dip 2 ) (5.2) 

To the non-cosmological approximation we neglect At against t and so, by (4.4), 

R 2 (t)/R 2 (t ) = t 2 /t 2 Q = 1 + 2At/t Q = 1 (5.3) 

the space-time becomes flat, the background expansion being insignificant. The non- 
cosmological approximation can simply be expressed as t$ 1 — > 0, in the sense that terms 
involving to in the denominator can be neglected. To this approximation G (hence 0) as- 
sumes its quasistatic value, for example 

G(t)/G(t ) = (t/t )~ n = 1 - nAt/to = 1 (5.4) 
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This is consistent with = O(t 2 ) -> by (4.13) and Gp m = O(t 2 ) -> by (4.11), i.e. 
the background becomes the vacuum and cosmological terms disappear. 

Now consider the sun as a local concentration of matter with density p(r) much higher 
than the cosmic background.. To the non-cosmo logical approximation we can solve the 
system of gravitational field equations (3.16) and (3.17) for g^ and <fi outside the sun by 
the quasistatic (t/t = 1 + At /to = 1) exterior solution 

Exterior (f)(t,r) = (f)(t ) = const (5.5) R^ v = (5.6) 

with the boundary condition at spatial infinity g^ v tending to, according to (5.3), i]^ u , that 
of flat space-time. Thus we see, as far as exterior solution are concerned, our theory reduces 
to Einstein's theory in the non-cosmological approximation. So our theory passes the crutial 
test about the advance of perihelion. 

VI. ELECTROMAGNETISM FOR VARYING G 

In this section we shall go to find the modification for varying G of the action integral 
for the electromagnetic field and the additional action integral for charged matter. We start 
from the large number equation (1.2). According to our approach of replacing evolving law 
by cosmic combined conservation law, for example see (2.4) (2.5) for I m for matter, which 
we consider it applicable also to the sun, the planet, even to a mass point, so 

M = G 1+1/n M, m p = G 1+1/n m p , m e = G l+1/n m e (6.1) 

and the left hand side of (1.2) is a fixed constant. Assuming c not evolving, we conclude 
that e 2 /m e oc t oc G^ 1 /™, so the cosmic combined charge 

e = G l/2+1/n e (6.2) 

conserves in the long. 

In Einstein's theory with constant G, the action integral for the electromagnetic field is 
(we follow again closely Dirac's book on general theory of relativity) 

I ern = -(167T)- 1 J F >lu F^^{-g)d i x where F^ = - (E6.3a,b) 

and the additional action integral for charged matter is (cf Dirac (29.4) (29.1)) 

I q = - J K^Td 4 x or = -e J K^dx^ (E6.4a,b) 
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for continuous distribution of charges or for a point charge. Recalling the rule of our modi- 
fication mentioned before in connection with I m and I m , i.e. mutiplying / with Q 1+1 / n and 
keep an eye on the conservation law, we get naturally, starting from (6.2) and (E6.4b), then 
(E6.4a) or (E6.3b), and finally (E6.3a), the following modified action integrals 

Ig = -e J H^dx^ with = G 1/2 k^ (6.4b,c) 
I q = - J njd^x withf = G 1/2+1/n f (6.4a,d) 

Ffiu ^n;u ^u;ti ^vifj, (6.3b) 

I em = -(lGn)- 1 J <P% u F^^/{-g)d A x (6.3a) 
We give the variations of (6.3a) 

SI em = y"[(-l/2)^% I/ + (47r)- 1 (0 2 Fn ; ^ 
+<$>emHW{-9)d 4 x (6.5a) 

with 

E; = 2 [-(4vr)- 1 F M(T F- + (IBtt)- 1 ^^^ (6.5b) 

$em = -(87r)-V^^ (6.5c) 

The variation of (6.4a) is 

5I g = - J [T^ + K^sf]d 4 x (6.6) 

In the comprehensive action principle 5I tot = 0, I to t — I g + Im + ^em + ^ we obtain from the 
coefficients of Sk^ the variational equation 

f = (47r)-V^W(-<?) = (4vr)- 1 [0 2 ^V(-^)] ) , (6-7) 
Hence the identity 

f in = 0, so J j dx x dx 2 dx 3 conserves (6.8) 

which reduces by (6.4d) to (6.2) for a point charge. In compliance with this conservation 
law, we take in (6.6) as familiar in (3.11) 
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After integrating by parts we obtain for th second part of (6.6) 

«m0 V i -fb v ), v d i x = J F, u ] b»d A x (6.10) 

For the elements of charged matter the equations of motion obtained from the coefficients 
of 6 M , by a change of notation similar to (3.9) (3. 10) 

1 = PeuW(-9) , Pe = G 1/2+1/n p e (6.11a,b) 
contains a Lorentz force term 

PmU^ + F^p e u v = (6.12) 

Dirac, in his book, uses a instead of our p e for the charge density. The variational equations 
for g^ v and <fi are 

N£ = -8ttT; - 8nE£, $ = -87r$ em (6.13a,b) 

The effect of the factor 2 in (6.3a) and (6.7) is familiar from Maxwell's theory; as seen 
from (6.7) in the case of Galelian metric 2 plays the role of relative dielectric constant 
and (0 2 )" 1 plays the role of relative magnetic susceptibility p^. Hence e^p^ = 1 which is 
consistent with our assumption that c does not evolve. 

VII. GEOMETRIC OPTICS 

We shall study in this section the propagation of light or electromagnetic wave in space 
time where the metric and <fi are given. As is well known in Einstein's theory of general 
relativity that light moves along a null geodesic, and this can be regarded as a good approxi- 
mation of geometric optics to wave optics for the electromagnetic equations in the Riemanian 
space-time. In this section we shall show that this is also true for the theory with varying 
G, following closely the smilar treatment 5 of the present author for constant G. 

From (6.3b) where the covariant derivative can be replaced by ordinary partial derivatives, 
we can eliminate the cosmic combined potentials and obtain 

Ffj, u ,\ +F u \, fl +Fx liu = (7.1) 
From (6.7) where charge-current is absent we have 

[<^F^y/(-g)], v = [^Wf-Sli- (7.2) 
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Among the four equations in the set (7.2) there is one differential identity namely the 
ordinary divergence of its left hand side vanishes identically. The same is true for the four 
equations in the set (7.1) which can be written in a form similar to (7.2) by introducing the 
dual of F^y defined by 

~41 _ ^23 _ 

F =F 23 ;F =F 41 ; 1,2, 3 cyclic (7.3) 

Hence the set (7.1), (7.2) contains six independent equations linear in the six dependent 
variables For light or electromagnetic waves, we need to consider solutions of the form 
of waves 

F^^f^sinS (7.4) 

where the phase S varies quickly in space or time at the scale of a wavelength or a period of 
light, but the amplitudes vary little at such scale, because the multipliers in (7.2) only 
vary at a scale much much larger. In contrast to the fast variable S we call / M „, (f),g pa all 
slow variables. Neglecting the derivatives of the slow variables in comparison with those of 
the fast variable, we obtain from (7.1) and (7.2) after removing the common factor cosS* 
and writing S, x = s x the equations 

f^sx + fu\s P + f^s v = (7.5) 

[^VV^/aM-s)]^ = i.e. s^hp = (7.6) 
where s 13 = g u/3 s u . In (7.6) the second equation follows from the first by contracting the 
latter with and removing the common factors 2 and g). These algebraic equations 
being homogeneous and linear in the amplitudes f pv , the condition for non-trivial solution is 
the vanishing of the determinant s x s x = 0. More simply this can be obtained by contracting 
(7.5) with s x and using (7.6) to obtain f fiu s\s x = 0. So at any point where there is light we 
must have 

sxs x = g Xfl s x s^ = (7.7) 

Let a real displacement along the light path be denoted by dx^, alomg the light path (7.7) 
always hold, so 

d(g^sxs^) = s x s^,„ dx v + 2s v ds„ = (7.8) 

Let a virtual displacement on the surface of constant phase i.e. S = const, be denoted by 
5x u . So we have 

= SS = S, u 5x" = s u 5x" (7.9) 
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Since these two displacement vectors are orthogonal to each other, we have also 

g^dafSa? = (7.10) 

On comparing the last two equations we see that the coefficients of 8x v there must be 
proportional. As is well known, we can choose the parameter p for the path of light x M = 
x^{p) such that (7.10) agrees with (7.9) and (7.8) is satisfied by the following two set of first 
order differential equations for the light paths or rays. 

g^dx^/dp) = s„, ds u /dp = -(l/2)g x > 1 , u sxs fl (7.11a,b) 

In Riemanian geometry, (7.11a,b) is the standard form for the differential equations of null 
geodesies, the two sets of fiest order equations being completely equivalent to the following 
one set of second order equations by eliminating s u with the help of (7.11a), 

d 2 x fl /dp 2 + T^dx" / dp){dx a / dp) = (7.12) 

and the equation (7.7) being equivalent to the null condition 

= g^s^ = g^g fip g U(T (dx f, /dp)(dx' 7 /dp) = ds 2 /dp 2 (7.13) 

As we have shown that light path in our theory as in Einstein's theory is the null geodesic 
and the metric outside the sun to the non-cosmological approximation is the same as that 
in Einstein's theory, so the crucial test about the deflection of light by the sun is also passed 
by our theory as good as by Einstein's theory. 

VIII. QUANTUM MECHANICS 

The fundamental natural constant in quantum mechanics is the Planck constant ft. In 
order to decide how it evolves, we need some action integral that contains it. It seems most 
reliable to consider the action integral for Dirac's relativistic theory of electron, because the 
generalization of Dirac's equation to curved space-time i.e. to general relativity has been 
considered by many early authors. Here we follow Schrodinger's treatment as summarized 
in my earlier paper p with minor changes in the notation (the electromagnetic potentials 
Ay, there is here written as and the companion there of the wave function ip is here 
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written as x )■ We consider hydrogen like atoms for simplicity. With constant G we have 

Id = J{(l/2)xG%iHd/dx^ + e K M 

+(l/2)[(-iHd/dx» + en,)x]G^ - X m^W(-g)d 4 x (8.1) 

where the four general gamma 4x4 matrices G M depend only on g^ u satisfying 

G^G" + G v G lx — 2g txv (8.2) 

and x an d 4> are respectively 1x4 and 4x1 matrices so that the chain product xG^"4> is 
a number. Here I D for the hydrogen atom includes the part previously denoted by I q . We 
have, cf (6.4a), for Dirac's electron 

f = -e X G"W(-<?) (8-3) 

the ordinary divergence of which vanishes as a consequence of the variational equations 
with respect to x an d ip, i.e. the Dirac equation 5Id/Sx — and its companion 5 Id /Sip = 
in curved space-time. 

U = X 5 j^-^=ixG^(-g)U (8-4) 

The variation of Id with respect to g a p can be obtained from the dependence of G M on g a p 
by using the formula derived from (8.2) 

dG»/dg aP = -(1/4)(GV" + GV M ) (8.5) 

According to our rule the modified action integral for varying G is 

Id = J{(l/2) X GV[(ihd/dx* + eKM 

+(l/2)[(-iHd/dx tM + eK fJ ,)x]G tM ijj - xm e ip}^/(-g)d 4 x (8.6) 

This confirms our previous relations (6.1) for m e and (6.2) (6.4c) for e k^. It shows further 
that the cosmic combined Planck constant 

h = G 1+1 ' n h (8.7) 

conserves in the long. 

Meanwhile we call attention to the expression of cosmic combined Coulomb potential 
k 4 due to a point charge, say Ze which, in the case of Galelian metric, is not Ze/ r but is 
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«4 = (0 2 ) 1 (Ze/r) = Ze/ (e^r) by solving (6.7) with // = 4 . Using (6.2) and (3.1), we verify 
that the Coulomb potential, only correctly written this way, evolves like (6.4c) 

k 4 = G 1/2 k a = G 1/2 Ze/r 
so ££ 4 = G 1+1/n en A = G 1+1/n [Ze 2 /r] (8.8) 

Hence the fine structure constant defined by (here c being restored) 

a = eVfofc) = e 2 /(hc) = a (8.9) 

remains unchanged like c during the evolution of the universe. 

We are interested in the atomic spectra (of hydrogen say) for light coming from extra- 
galactic nebulae or from the sun. To obtain the energy levels by quantum mechanics we 
shall follow closely the treatment in my earlier paper |5j in Einstein's theory. We use the 
local Galelian metric obtained by a coordinate transformation, which can always be done if 
we treat the metric g^ v as constants, assuming their value at the position of the atom at the 
time of emission or absorption of the light quantum. For simplicity consider the diagonal 
case 

9u = -(A*) 2 , % = 1, 2, 3, g u = +(A 4 ) 2 , = i/) (8.10) 

Then we have 

G l = (\)-y,i = 1,2,3, G 4 = (A 4 )-y (8.11) 
where 7 M are the 4x4 gamma matrices familiar for special relativity, defined by 

yy + yy* = 2^ (8.12) 

Then Dirac's equation as obtained from (8.6) by varying \ becomes 

{7 4 (A 4 )~ 1 [^(9/9t) +S£ 4 ] +j j (X j )- 1 m(d/dx j ) - m e }tfj = (8.13) 

where the vector potential vanishes and the scalar potential due to the nucleus at rest at 
the origin is to be obtained by solving (6.7) with /i = 4, namely 

[0 2 F 4 V(-S)]„= 4tt(Z?)5(x 1 )(5(x 2 )5(x 3 ) (8.14) 

In the region of atomic scale </> may also be treated as constant, so the left hand side of 
this equation becomes simply ^g^g-^A^AsA^K/i),^ summed over j. It is convenient to 
introduce the scaled variables separately 

y 3 = XjX 3 so 5(y 3 ) = 5(x 3 )/\j j = 1,2,3 not summed (8.15) 
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Then (8.14) becomes in the scaled variables 



Vj(£ 4 ) = -H^^WWW) (8.16) 

the solution being proportional to the inverse radial distance r y of the scaled y's 

k 4 = X^Ze/^r- 1 (8.17) 

With the help of (8.17), the equation (8.13) can be written in the scaled variables 

{^[{X^-HHid/dt) + ZeV -2 ^ 1 ] + j j [ih(d/dy j )} - m}^j = (8.18) 

or after removing the common factor Q 1+1 / n (cf.(8.3) (6.1) (6.2) (3.1)) 

{^[{X 4 )-Hh{d/dt) + (Ze 2 )r y l + 7 j [iH(d/dy j )} - m}^ = (8.19) 

The energy value including the rest energy defined by 

ih(d/dt)ip = Eip (8.20) iH(d / dt)ip = Eif; (8.21) 

is easily obtained by comparing with that of Dirac's equation in his book on quantum 
mechanics.. Using the quantum numbers n r and j we obtain (with c restored and using 
(8.5) for A 4 ) 

E nr „ 3 /mc 2 = E nr Jmc 2 = (g 44 ) 1/2 {l + ff ^ V 1 ' 2 (8.22) 

For a transition between two quantum states A and B with quantum numbers n r (A),j(A) 
and n r (B), j(B) the frequency u A ^ B of the atomic spectral line is 

va-+b = (E A - E B )/h = (E A - E B )/h = (g 4 4) 1/2 (vA^B) Q M (8.23) 
where the ordinary quantum mechanical value of the frequency is given by 

(^s)qm = ^-{1 + - ff 72 V 112 \i (8.24) 
27rA c (n r + y/f - Z 2 a 2 ) 2 

Since both the Compton wave length X c = H/mc = h/fhc = X c and the fine structure 
constant a = a remain unchanged, we see from (8.24) that the frequency of the atomic 
spectral line does not change during the evolution. Also the Rydberg constant -Roo = 
mc 2 a 2 / (4irh) = remains unchanged. So does the mass ratio e.g. m p /m e = m p /m e . 
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The relation (8.23), it seems, can be generalized to all energy levels of complex atoms and 
molecules. 

We have shown in (8.22) that the gravitational effect on the frequency of the spectral 
line during emission or absorption in our theory is given by the factor (g^) 1 ^ 2 , which is the 
same in our theory as in Einstein's theory to the non-cosmological approximation. So the 
crutial test about the gravitational red shift of the sun is also passed by our theory as by 
Einstein's theory. 

We note in the Robertson- Walker metric using cosmic time with g 44 = 1 the frequency 
of the spectral lines emitted anywhere at any time is by (8.23) always the same. 

IX. HUBBLE RELATION 

In this section we shall investigate the relation between the cosmic red shift of the spectral 
line received from a cosmic distant source and the distance of that source from an observer 
here at present. We use the Robert son- Walker metric (5.2), so the frequency of the spectral 
lines for the distant and nearby object are the same by (8.23), namely uqm just when they 
are emitted. But there is a change of the frequency of the spectral line of the distant source 
during its propagating towards the receiver according to (7.11b). 

We denote quantities referring to the distant source by the suffix s and those referring 
to the receiver or observer here at present by the suffix 0. For simplicity we use geometric 
optics and radial rays, putting d6 = dip = sg = = 0, the phase of the spherical waves 
being, as usual, S = 2n(ut — r/A) = 2itv{t — r/f$), with s 4 = 2ttu, s r = —2-k/X. The null 
condition gives the velocity of light (3 for the metric (5.2) 

Xv = P = R(t )/R(t) (9.1) 

So at the time t s the wavelength of the spectral line emitted by the source is 

\,VQM = R(t )/R(t s ) (t s <t ) (9.2) 

For the metric (5.2) , being homogeneous, the wave length of this spectra line keeps constant 
by (7.11b) during the propagation, till the line being recieved and compared at time to with 
the same line emitted near by with 

Ao^qm = 1 (9.3) 
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The cosmic redshift z is defined by z = (A s — A )/A , so dividing (9.2) by (9.3) we have 

l + z = \ s /\ = R(t )/R{t s ) (9.4) 
The optical distance A the light propagated is obtained by integrating dr from t s to t 



r*° R(t 



A = J. dr = I, m dt (9 5) 

If this is much larger compared to the linear dimension of the source the inverse square law 
holds for the intensity and can be used to measure the distance. 

For matter dominated cosmological model we have R(t )/R(t) = t /t, so from (9.5) and 
(9.4) we obtain 

A = t log(* /*«) = *o log(l + z) (9.6) 
This becomes the Hubble relation 

z = H A with H Q = (to)- 1 for z « 1 (9.7) 

but for larger z we have from (9.6) 

z = exp(# A) - 1 (9.8) 

X. ELECTROMAGNETIC RADIATION 

In the next section we shall supplement the matter dominated universe with electromag- 
netic radiation. Owing to the appearance of the factor 2 in the action (6.3a) we need study 
the electromagnetic waves a bit more in order to obtain the expressions corresponding to 
(6.5b) and (6.5c) for electromagnetic radiation. 

We write at first the electromagnetic wave equations in a form familiar from Maxwell's 
theory. Free from source (6.7) splits into 

divD = 0, (lO.l)D = [0V(-^)](^ 41 ^ 42 ^ 43 ) (10.1a) 
curlH -dD/dt = (10.2) H = [0V(-#)](F 23 , F 31 , F 12 ) (10.2a) 

The other pair of Maxwell's equations come from (6.3b), from which we obtain 
+F\ ll , v +F uX ,ix = which splits into 

divB = (10.3) B =(F23,F 3 i,Fi 2 ) (10.3a) 
curlE+dB/dt = (10.4) ~E = (F u , F 24 , F u ) (10.4a) 
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Comparing (6.5c) with (6.3a), the integrand of the latter being equal to two times (B ■ H — 
E ■ D), we obtain 

$ em = -(B ■ H - E ■ D)/[^{-g)] (10.5) 

For the Robertson- Walker metric with k — 0, it is convenient to use Cartesian coordinates 
to write 

ds 2 = dt 2 - f 2 (dx 2 + dy 2 + dz 2 ) with / = R(t)/R(t ) (10.6) 
Then we have from (10.1a) to (10.4a) the slowly varying permittences e and /i 

D = eE, (e = [0 2 / 3 ]/- 2 ) (10.7a) 
H = ^B, (/.- 1 = [0 2 / 3 ]/" 4 ) (10.7b) 

We treat the electromagnetic radiation as superposition of plane electromagnetic waves. 
For a plane electromagnetic wave propagating along the direction denoted by the unit vector 
s , the only fast variable is the phase S = ~r • s — j3t where (5 = 1/y/eJI, like e and fi, 
vary slowly in time and will be treated as constants. Then (we follow here Born and Wolf's 
treatment 6 ) (10.2) and (10.4) becomes 

~E = -^Je(~s xH),H = ^e[ii(~s x E) (10.8) 

Hence E , if, and "s > form a right-handed orthogonal triad, and we have 

eE 2 = fiH 2 orD-E-H-~B=0 (10.9) 

So by (10.5) we have $ em = for a plane electromagnetic wave, and by superposition of 
plane waves we obtain for electromagnetic radiation 

$em = (electromagnetic radiation) (10.10) 

For electromagnetic field where charge-current density vanishes we have the differential 
identity, 

E^ + $ em ^ = (10.11) 

which can easily be verified from (6.5b) and (6.5c). Where charge-current density does not 
vanish, this identity should be replaced by 

f^ + E^ + ^.,^0 (10.12) 
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so that the terms involving Lorentz's force from T^. v and E v cancel. For electromagnetic 
radiation (10,11) simplifies by (10.10) to 

E v = (electromagnetic radiation) (10.13) 

For homogeneous and isotropic electromagnetic radiation the superposition of plane waves 
propagating in all directions leaves only the diagonal elements of E v non-vanishing Since 
from (6.5b) we have zero trace before the superposition, we conclude that 

E\ = El = El = = "Pr/3, E\ = p r (10.14) 

It is well known that with (10.14) and the metric (10.6) that (10.13) can be integrated to 
obtain 

Pr = Pr(t) = Pr(to)R 4 (t )/R 4 (10.15) 

XI. MATTER PLUS ELECTROMAGNETIC RADIATION 

We now consider homogeneous cosmology including matter and electromagnetic radia- 
tion, I tot = I g + I m + I em and write according to our theory with varying G the equations 
for the comprehensive variational principle 5I tot = 0. 

The variational equation for remains unchanged as (3.17) because of (10.10). The 
variational equations for g^ u are 

Nt = -8n(p rn + p r ), Nl = 8vrp r = 8np r /3 (11.1) 

We give these equations in full, cf. (4.8), (4.9) and (4.10) 

. 2 

$ = 6[R/R + R /R 2 }(j) 
+ 8[0 + 3(i?/ J R)0] = (11.2) 

Nl = (f) 2 [-R 2 /R 2 - 2R/R] - (<f> 2 )" - 2(<f) 2 yR/R 

. 2 

+40 = (87r/3)p r (t )R 4 (t )/R 4 (11.3) 

.2 .2 

Nl = <j) 2 [-3R /R 2 ] - 3 (0 2 )' R/R- 40 

= -87r\fl m (t )R 3 (t )/R 3 + p r (to)R 4 (t )/R 4 } (11.4) 
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The differential identity, cf. (3.6), + $0 ;4 = shows that among the three equations 
(11.2,3,4) only two are independent, the condition for compatiability on the right hand side 
being separately guaranteed for matter and for electromagnetic radiation. The combination 
N" + $0 cf (3.18) remains unchanged and can be integrated to give, cf. (4.2), 

(0 2 )- J R 3 / J R 3 (t o ) = -8np m (to)t (11.5) 

which can be used in place of any equation among (11.2,3,4). We may take as the two 
independent equations (11.4) and (11.5), and use the new variables 

r = t/t ,S = R/R(to),r, = 2 /0 2 (*o) (11.6) 

and parameters 

87T^p m (t o )/0 2 (to) = tort 2 G(t )p m (t ) = A (11.7) 

Pr(to)/Pm(to) = Pr(t )/ Pm (t ) = E r/m ~ (1(T 2 to 1(T 4 ) (11.8) 

Then (11.4) and (11.5) become, with circle denoting d/dr, 

Svf/e - 3^|/<e - V/V = -A/e - e r/m A/? (11.9) 
^ 3 = -At (11.10) 

The solution for matter dominated universe given in §4 is for e r / m = that £ = r, rj — 1/r, 
and A — 1. This solution satisfies not only the "initial" conditions at r = 1 that £ = 1 
and i] = 1 but also the "final" condition at r = that £ = 0. The extra condition serves 
to determine the parameter A. We shall leave the mathematical problem of solving these 
equations (11.9), (11.10) but only note that the asymptotic solution of these equations for 
large £ is that of matter dominated universe. Only with the help of the asymptotic solution 
which is simple enough to reveal the cosmical combined conservation laws that we find the 
way to modify or to generalize Einstein's theory of general relativity to work to cosmological 
time. For small values of r, at r ~ s r / m , the deviation from the asymptotic solution will 
be appreciable. For smaller values of r refinement of the action integral for matter with 
corresponding refinement of need to be considered to take account of the internal energy 
and pressure, which at thermal equilibrium can be treated according to the principle of 
statistical mechanics. 
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XII. STATISTICAL MECHANICS 



The fundamental natural constant in statistical mechanics is the Boltzmann constant k B , 
in this section we shall see how it evolves in cosmological long time. For this purpose we 
need to consider for example black body radiation where k B appears in Planck's formula 
and in the energy density of black body radiation. 

f°° 8nhu 3 /c 3 , 7T 2 fc* . ^ 
Pr= — 77-717^ dv=—-§T* (12.1) 



Iq exp(hu/k B T) — 1 15c 3 h 3 
Comparing this with the cosmic combined expression 



^hu 3 /c 3 , 7T 2 (k B f 



<,7l his c 



p r = / dv= -^-^T 4 (12.2) 

Jo exp{hu/k B T) - 1 15c 3 (Hf 

and noting that like p m (cf. (11.1) and (2.5)) we must have p r = G 1+1 ^ n p r , as we know that 
electromagnetic radiation can be regarded as matter, namely photons. From this we obtain 
from (12.1) (12.2) and (8.7) the similar relation 

k B = G 1+1/n k B (12.3) 

Thus we have the general rule for c.c. energy, cf. (2.5), (8.7), (12.3), and also cf.(8.8) for 
the c.c. Coulomb energy, 

E = G 1+1/n E, (E = mc 2 M,k B T,ZR 4 ) (12.4) 

where c 2 ,u,T do not change. Since energy, momentum and stress form a four dimensional 
tensor, the above relation (12.4) holds too for the momentum p x ,P y ,Pz and the pressure P, 
while velocity u x ,u y ,u z and volume V do not change. This makes the Boltzmann factor 

exp(-E/k B T) = exp(-E/k B T) (12.5) 

the same and the counting of phase cells 

Vdp x dp y dp z / (h) 3 = V dp x dp y dp z /h 3 (12.6) 

also the same, so for canonical ensemble of systems at local thermal equilibrium, the partition 
function (Zustandsumme) Z = Z and hence the distribution function is the same with 
or without the cosmic combination factor. For example, for Bose-Einstein statistics, the 
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distribution function of photons is given by, cf (12.1) and (12.2), where the factor two 

accounts for the two polarizations 

(2)47rz/ 2 /c 3 (2)47rz/ 2 /c 3 



exp(hv/k B T) -1 exp(hv/k B T) -1 



;i2.7) 



We shall not go into the details of the early universe but suffice it to note that at high 
temperature the internal energy and pressure must be taken into account. In Einstein's 
theory the action Integra! /. (E2.3a) for matter withont internal stress is replaced by that of 
perfect fluid Ifi, {we follow here Fock's treatment in his book pj space time and gravitation, 
eq. (48.28) et seq.}, the cosmic combined expression in our theory is given accordingly by 

= - J p>M + t m )V(-9)d 4 x (12.8) 

Here e m denotes the c.c. internal energy per unit c.c. mass p m of the ideal fluid and is 
considered as a function of the latter only, as we are dealing with adiabatic processes in 
which the c.c. entropy per unit c.c. mass, denoted here by s m is kept constant, i.e. 

de m + Pd(l/p m ) = TcB m = (12.9) 

To find the variation of Ifi we express p m in terms of by (2.4a) and obtain 

51 f i = - [ -j^-5(PmV ::: g)d A x 
J dp m 

f dF 
+ {Pm-^ - F){5^-g)d A x 

J dp m 

= -(1/2) J (fpg,„W(-g)d*x 
I dF 

- / —ujp^d^x (12.10) 

J dp m 

with 

~ dF dF 

T% = Pm^u^ - (p m ^— - F)g^ 

1 dp m dp m 

= [p m {l + t m )+P}u^u u -Pf v (12.11) 

Using (3.11) in compliance with the modified constraint (2.3b) we obtain from the second 
part of (12.10) the following equations of motion for the ideal fluid 

_ dF „ dF ~ ,dF 

Pm-T^U U K1/ + p m U U^—^j-u - p m (-p-) ; ^ 
dp m (LPm dPrn 

= \p m (l+e m )+P]u v u KV + u"u li P v ,-P ;ii = (12.12) 
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One can verify the covariant divergence identity 

~ ^ dF ^ dF _ dF 

T fi.w = (pmU u ) ]U -^u^ + p m u v {-—u il ), v - (/0 m -p- - F) ; ^ = (12.13) 

rfPm CtPm dp m 

with the help of (12.12) and (2.3b) (2.4b). We note that Ifi, like I m , does not involve 0, 
hence contributes nothing to the variational equation with respect to <fi. 

In the specific internal energy e m , being energy per unit mass, according to (12.4), the 
cosmic combination factor will cancel. Thus we expect e m = e m . This can be easily verified, 
e.g. for an ideal gas at low temperature. With Maxwell's law for the distribution of velocity 
exp[— m(ul+Uy+ul) / (2kBT)}du x du y du z in which rh and k B can be simultaneously be replced 
by m and k B , the translational internal energy is p m e m = 3n m k B T/2 while p m = n m rh 
where rh being the c.c. conservered mass of one molecule, and n m the number density of 
the molecules. Similarly ? m = e m will hold for relativistic gas at high temperatures. 



XIII. SUMMARY AND REMARK 



From our theory with varying G we have sysmatically found the cosmic combined natural 
constants rh, H, k B ,e 

m/m = H/H = k B /k B = [G/G ] 1+1/n , e/e = [G/G ] 1/2+1/n (13.1) 

which remain to be constant in the long. Our present fundamental principles of quantum 
mechanics and statistical mechanics work in the long with rh, H and ks- It is due to the 
old age of the present universe that the combination factors (7 1+1 / n involved in these c.c. 
constants vary very slowly by now that we take m, H and k B as natural constants. The law 
of electromagnetism in the long differs from the law at present by a dielectric and magnetic 
permittivity e = p^ 1 = [G/Gq] 1 ^. The velocity of light c, the Compton wave length for 
a particle H/fhc, the fine structure constant ^/(eHc) and the proton-electron mass ratio 
rh p /rh e all remain unchanged during the evolution. 

For phenomena occurring outside the gravitating body like the crucial tests our theory 
to the non-cosmological approximation reduces to Einstein's theory. Difference arises for 
phenomena occurring inside the gravitating body, including the case of cosmology. E.g. for 
matter-dominated universe our theory with varying G yields for the expansion law R oc t 
between the radius and the age of universe so as to be in agreement with the empirical 
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large number equality GM oc t, while in Einstein's theory with constant G and conservation 
of mass we have R oc t 2 / 3 with GM = constant of course. In our theory with varying 
G , a tensor term automatically arises from the spatial and temporal derivatives of G. 
This tensor automatically vanishes in a region where G is constant, as is the case with the 
exterior solution of the sun when we consider the crucial tests of general relativity. In the 
problem of non-static homogeneous cosmology this tensor with non-vanishing components 
A\ = A 2 = A 3 = l/t 2 , Aj = 2/t 2 (Or less probably A\ = A 2 = A| = 2/t 2 , A\ = l/t 2 ) 
takes the place of cosmological constant. It seems desirable to make an analysis of the 
observational data from the beginning according to the present theory with varying G. It 
would be also interesting to apply the theory with varying G to the interior solution of some 
high density objects where spatial variation of G is appreciable. 
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